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1 Introduction 



We illustrate a formalism that allows to match the next-to-leading order cross 
section for the photoproduction of heavy quarks to the cross section obtained 
by resumming logarithms of p T /m to the next-to-leading accuracy, thus giving 
a sensible prediction for any value of Pt- We present a comparison between our 
predictions and HI and ZEUS data. 
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The computations of the cross sections for the production of open heavy 
flavours at fixed-order (FO) in perturbative QCD have been so far carried 
out at the next-to-leading order (NLO) accuracy for all possible production 
<*i ' mechanisms. For single-inclusive observables, when the transverse momentum 

ft, p T of the quark is much larger than its mass m, terms containing log(p T /m) 

O ' become large in the FO cross section, which is thus no longer reliable. In 

1} , such a case, techniques are used that allow the resummation to all orders of 

i-G ■ the leading (LL) and next-to-leading (NLL) towers of logarithms of p T /m; in 

order to achieve the resummation, the heavy quark is effectively treated as 
massless. 

r% | The resummed (RS) and FO cross sections are therefore relevant to corn- 

ed ' plcmcntary regions in p T . The problem then arises whether to compare the 

data to FO or to RS predictions, since the inequality p T 3> rn cannot be 
unambiguously translated into a quantitative statement. It is thus desirable 
to write the single-inclusive cross section in a form that is sensible in the 
whole p T range, that is, which interpolates between the FO result, relevant 
to the small- and intermediate-p T regions., and the RS result, relevant to the 
large-p T region. This is the aim of ref. Ill and ref. □, relevant to hadro- and 
photoproduction respectively. In this letter, we shall concentrate on the lat- 
ter, reviewing the formalism and presenting phenomenological predictions for 



charm production at HERA. 
Our master formula reads: 

FONLL = FO + (RS - FOMO) x G(m,p T ) , (1) 

where FONLL (for Fixed Order plus Next-to-Leading Logarithms) will give 
sensible predictions in the whole p T range, and FOMO is obtained from FO by 
letting to zero all the terms suppressed by powers of m/p T . The subtraction 
of FOMO from RS in eq. ([j]) is necessary to avoid double counting, since some 
of the logarithms appearing in RS are already present in FO. To be more 
precise, FONLL will have the following features: 

• All terms of order a cm a s and a om a 2 are included exactly, including mass 
effects; 

• All terms of order a cm a s (a s logp T /m) 1 and a em aj, (a s log p T /m) 1 are 
included, with the possible exception of terms that are suppressed by 
powers of m/p T . 

Finally, the function G(m,p T ) is rather arbitrary, except that it must be a 
smooth function, and that it must approach one when m/p T — > 0, up to 
terms suppressed by powers of m/p T . In what follows, we shall use 

G(m,p T ) = 2 f\ 2 , (2) 

with c a free parameter. The practical implementation of eq. (ffl) is rather 
involved, especially in the photoproduction case; all the details can be found 
in ref. B 



2 Charm production at HERA 

In this section, we shall compare the predictions based upon eq. (|l|) to the 
experimental data for D* meson photoproduction obtained at the HERA col- 
lider by the HI U and the ZEUS □ collaborations. We set m=1.5 GeV, and 
the renormalization and factorization scales equal to the transverse mass of 
the quark, yp% + m 2 . The parton densities in the proton (photon) are given 
by the CTEQ5M1 (AFG) set. We set A^ = 226 MeV, as constrained by 
the CTEQ5M1 set; this value is almost identical to the central value of the 
PDG global fit. Finally, in order to simulate the hadronization of the bare 
c quarks into c-flavoured mesons, we use the Peterson function, and multi- 
ply the results by the probability of a c quark fragmenting into a D* meson, 
Pc^d* ■ In order to match the experimental conditions, we also need to apply 
cuts on the variables y and Q 2 , that enter the Weizsacker- Williams function 



(the latter variable is the upper limit of the photon virtuality squared). In 
particular, in the case of the ZEUS data, we have: 

0.187 < y < 0.869, Q 2 = 1 GeV 2 . (3) 

HI has two data samples, corresponding to the two electron taggers (ETAG33 
and ETAG44) used in the analyses: 

0.29 < y < 0.62, Q 2 = 0.01 GeV 2 : ETAG33; (4) 

0.02 < y < 0.32, Q 2 = 0.009 GeV 2 : ETAG44. (5) 

We also point out that ZEUS use P C ^ D » = 22.2%, and HI use P C ^ D * = 27%. 

In figure [j|, ZEUS data are compared to our theoretical predictions. The 
solid line is obtained from eq. (|l|) , by setting c — 5 and the Peterson parameter 
e = 0.02. The dotted (dot-dashed) line is obtained with c = and e = 0.02 
(c = 5 and e = 0.036). Finally, the dashed line is the FO result (at NLO 
accuracy), with e = 0.02. It appears that none of the theoretical curves 
describes the experimental data particularly well, both in normalization and 
in shape. The impact of unknown power-suppressed terms, and of terms 
beyond NLO in the perturbative expansion of the resummed cross sections, 
can be roughly estimated by looking at the difference between the matched 
results obtained with c — 5 and c = (solid and dotted curves respectively). 
This difference is non negligible, but it does not help in understanding the 
discrepancy with the data. On the other hand, the comparison between theory 
and HI data, fig. |[ appears to be satisfactory. The same pattern is found 
in the rapidity and pseudorapidity spectra: HI data are well reproduced by 
QCD, ZEUS data are not. A more detailed comparison between theory and 
HERA data will be presented in a forthcoming publication, where a complete 
study of the dependence upon the parameters entering the calculation (such 
as the mass and the scales) will be carried out. 

One of the authors (S.F.) wishes to thank C. Grab and L. Gladilin for 
useful discussions on HI and ZEUS experiments, and CERN TH Division for 
the hospitality during the course of this work. 
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O: ZEUS data 
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dashed: FO, e = 0.02 
solid: FONLL, c=5, e = 0.02 
dotted: FONLL, c = 0, e = 0.02 
dotdash: FONLL, c = 5, e = 0.036 
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Figure 1. p T spectrum of D* mesons, in the visible region of the ZEUS experiment. ZEUS 
data are compared to the theoretical predictions. See the text for details. 
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Figure 2. p T spectrum of D* mesons, in the visible regions of the HI experiment: ETAG44 
(left panel), and ETAG33 (right panel). HI data are compared to the theoretical predictions. 
The patterns of the theoretical predictions are as in fig. hi 



